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The baryon structure is investigated in a covariant diquark-quark model. In this 
approach baryons emerge as relativistic bound states of a constituent quark and a 
0"*" or 1+ diquark. After solving the Bethe-Salpeter Equation for the scalar diquark 
quark system in ladder approximation we couple various external currents to the 
constituents of the baryon to probe its internal structure. The quark and the 
diquarks are assumed to be confined which is implemented by suitable choices for 
the propagators. This leads to nontrivial vertex functions between the constituents 
and the external current. Baryonic matrix elements are then evaluated to extract 
observable formfactors. 



1 Introduction 

The main goal of our work is to understand the baryon structure up to a mo- 
mentum transfer of « 1 — 2 GeV. This momentum regime is particulary 
interesting because of the interplay between hadronic degrees of freedom on 
the one side and their intrinsic quark structure on the other side. To describe 
hadrons as bound states of quarks nonperturbative methods of QCD are un- 
avoidable, whereas for large the perturbative results of QCD should be met. 
So we are looking for an interpolating model to describe baryon structure at 
the intermediate energy region. 

From the experimental side there is also an interest in this kind of cal- 
culations, because proposed and already started experiments (e.g. at COSY, 
ELSA, TJNAF) explore the baryon structure to a very high precision in this 
momentum regime. 

This workshop is devoted to diquarks as a tool to parametrize complicated 
or unknown structures. We will exploit this philosophy in the following. Since a 
fully relativistic Faddeev equation which determines the nucleon as bound state 
of three quarks (bound by some kind of gluon interaction) is almost impossible 
to solve it is worthwhile to study the approximation where two quarks are 
bound to a diquark and interact with the third quark through quark exchange. 
Such a picture of a baryon has been derived. by path integral techniques in the 
context of the Nambu-Jona-Lasinio model u and in the Global Color Model □. 
* Supported by COSY under contract 41315266. 
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2 Bethe-Salpeter Equation for Baryons 

The basic assumption of our model is that a baryon is a bound state of a 
constituent quark and a scalar (0^) or axialvector (1^) diquark which are 
bound due to quark exchange EJ. The corresponding Bethe-Salpeter equation 
(BSE) for the nucleon is then given by Q 

HP,P) = J ■0^^g{-k-p)^^S{^k~p)^^S{\P + k)D{\P-k)^{P,k) 

+ j -0^gi-k-p)rS{-k-p)j5SC2P + k)Df^'^C^P-k)^'^{P,k) 

^^{P,P) - j (045(-fc - P)l^oS{-k - p)rS{^P + k)D{\P - k)<i>iP, k) 

+ / ~ ~ Ph'^SilP + k)D^P{}P - k)^P{P, k) 

(1) 

where $(P, p) and ^{P,p)'^ are the baryon (scalar and axial vector) vertex 
functions with amputated quark and diquark legs. The quark propagator is 
denoted by S{q) and the diquark propagator by D{q) and D'^'^{q) respectively. 
g{q) describes a possible extension of the diquark quark vertices. A solution 
of the BethcpSalpeter equation which assumes a static quark exchange can he 
found in ref. Although 1+ diquarks are partially included in our calculationa, 
we will restrict ourselves in the following to scalar diquarks so that only the 
first line of equation (1) survives. 

Because $(P, p) describes a spin i particle it is useful to write it as a 
product of an amplitude '^{P,p), which depends on the_total and the relative 
momentum of the nucleon, and a Dirac spinor u{P, S): U 



$(P,p) :=*(P,p)w(P,5). (2) 

After inserting this ansatz in the Bethe-Salpeter Equation (^, we multiply 
with the adjoint spinor u(P, 5*) from the right hand side and sum over the 
nucleon spins. It is the inimidiately seen, that only the projections 

X(P,P):=*(P,^)A+, A+ = i^!^,I±^ (3) 



"We use a Euclidean space formulation with {7^1, 7iy} = 25^1/, = 7^, 
p ■ q = ~ Pi'ii assume isospin symmetry. 
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to positive energy survive. The projected amplitude x{PyP) is then by con- 
struction an eigenfunction of A+, 



X{P,P) - ^{P,p)A+ = x(P,p)A+. (4) 

This leads to the possibility to decompose further the Dirac structure of x ■ 
The decomposition 

X{P,P) = i(l - + U-n-PT - 2k^l-Pl-P- 1-P1-P))S2, (5) 

is the only one compatible with equation . p!^ — p^^ + denotes the 

transversal relative momentum. The still unknown functions Si{P^ ,pP,p'^) 
and S2{P^ ,pP,p^) are of course determined by the solution of the BSE. 

In the rest frame of the bound state (where the actual calculations are 
done) we obtain 

. / lSi{P,p) \ 

x.,.{P,P)-[ff^S2{P,p) j- 

It is now obvious that the first two columns of the 4x4 matrix correspond 
to spinors with positive energy, spin up ("large" component) and spin down 
("small" component) respectively. In the calculation of baryonic matrix ele- 
ments we will restrict ourselves, in a first approximation, therefore to Si. 

Finally we expand the scalar functions Sj (P, p),j = 1, 2 in terms of Gegen- 



bauer Polynomials Q 



oo 

S,{P,p) ^y^^^fi{^)Cl{ L,' ^ \ (7) 



Sr^ :k li I P ' P 

— 



and solve the Bethe-Salpeter Equation in the rest frame of the nucleon either by 
diagonalization or by iteration. In this way, the bound state mass is determined 
by the eigenvalue of the integral equation, while the expansion coefficients 
correspond to the eigenvectors. 



2. 1 Parametrization of confined constituents 

Since neither quarks nor diquarks have ever been observed as asymptotic states 
in experiments, there is not much justification necessary to include this feature 
of QCD in a calculation of baryon properties. 

We implement confinement by using propagators which have no poles. 
In case of the quark propagator this has to be understood as an effective 
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parametrization of the behavior obtained in Dyson-Schwinger studies of qcdI 
The quark self-energies which appear in 

S{p) = ~i-fpa,{p) + a,{p) = [^7•M(P) + B{P)]-^ (8) 

are then choosen as 

l-exp(-(l + ^) 

(J^{p) = 5— 5 2_ and asip) = mqay{p). (9) 

p'^ + 

Obviously the exponentials remove the mass poles however, thereby introduc- 
ing aiL essential singularity at = —00. 

In refO it has been shown that confined diquarks can be obtained when one uses 
an appropriate irreducible quark-quark kernel (beyond ladder approximation). 
We simulate this by using 

Dip) = ^ (10) 

where ^ 

Z(p) = l-exp(-(l-f ^)). (11) 

Note that an elementary scalar diquark would correspond to Z = 1. With such 
a prescription of confined constituents we especially get rid of the unphysical 
quark-diquark thresholds which plague baryon calculations within the pure 
hadronized NJL-modelEl. 

To estimate the effects of an intrinsic diquark structure we did the calculation 
with a pointlike quark-diquark coupling 

9(k)^go (12) 
as well as with a diquark formfactor (see also ref. 0) 

A2 

gW -g0 fc2^A2 ' ^^^^ 

where A is of the order of 2mg and which simulates the diquark structure in a 
crude way. 

2.2 Numerical results 

In Fig. 1 the numerical results for the dominant expansion functions fl without 
and with a diquark formfactor are shown. The parameters here and in the 
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Figure 1: Numerical Solution of the Bethe-Salpeter equation. The results in the left figure 
are obtained with a diquark formfactor, in the right figure with a pointlike coupling. 



following are fixed to rriq — = 0.5 GeV and the coupling is chosen to 
obtain Mb = 0.95 GeV. In both figures the functions fall off rapidly for large 
relative momenta, but the asymptotic behavior is quite different: Assuming 
a diquark structure leads to a narrower amplitude in momentum space. This 
will also have an effect for the observables, where a possible diquark structure 
just enters via the nucleon vertex function. We observe that in both cases the 
expansion in terms of Gegenbauer Polynomials converge rather rapidly. 

3 Hadronic matrix elements 

3. 1 Mandelstam 's formalism 

Matrix elements of operators in the Bethe-Salpeter approach are calculated 
according to Mandelstam's formalismllil 

=/(0j/ 0^^iP'^P)^oip'^P'-^P^PMP^P)- (14) 
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Note that P and p [P' and p') are the total and relative momenta of the 
diquark-quark states before and after the interaction; ip is the normalized 
diquark-quark Bethe-Salpeter wave function related to the vertexfunction by 
multiplying with the quark and diquark propagator. The 5-point function 
describes how the external current couples to all internal lines of the diquark- 
quark system. 

In the following we work in a generalized impulse approximation, where 
we consider only the coupling to the quark and to the diquark, 

T^{p\ P'-p, P) = Viip', P';p, P) + T%{p', P'-p, P), (15) 

and neglect the coupling to the exchanged quark. 

3.2 Vertexfunctions of the constituents 

In case of electromagnetic (e.m.) formfactors, we still have to know the quark- 
photon and the diquark-photon vertexfunctions. Using gauge symmetry, which 
is manifest in the Ward-Takahashi identity, 

+ = S^^iP + Q) - S^^ip) (16) 

one observes the connection of the longitudinal part of the vertexfunction and 
the inverse propagator. 

A vertexfunction which solves the above identity for the quark is the Ball-Chiu 
vertex t3 

-^{Bip^)-Bie))}, (17) 



and in case of the diquark 



Tf^ip,k) = (p + k),^^^^^-^^. (18) 



These vertexfunctions are totally determined by the self-energy functions of 
the corresponding propagators. 

4 Electromagnetic nucleon formfactors 

Within our approximations we decompose the matrix element into a quark and 
a diquark part 

MQ^) - J^(Q') + QdJ'(Q'), (19) 
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Figure 2: Electromagnetic Formfactors of the nucleon: The soUd Unes are the proton formfac- 
tors, the dotted hnes the neutron formfactors; the thin curves are the results with a diquark 
formfactor, whereas the thick curves are the results assuming a pointlike quark-diquark 
coupling. 



(20) 



with the charge matrices in isospin space 

The quark part is given by 

J',{Q') = J-0^m,Pf)Sq{k+)iT^''{k+,k.)S,{k.)S4ka)x{P^^ (21) 
whereas the diquark part is calculated with the loop diagram 

'^'^^x(^/,p/)5rf(fc+)ir^(fc+,fc_)5d(fc_)S,(fc,)x(Pi,Pi). (22) 



From Lorentz invariance and discrete symmetries it is also known that the 
onshell e.m. nucleon current can be decomposed into 



{N{Pf,Sf)\j^\NiP,,S,)) 

u{Pf,Sf)[lMB{Fe{Q^) - Fm{Q^))^+Fm{Q^)-f^]u{Pi,S,) 



(23) 



Note that in our formalism the Dirac spinors appearing on the right hand side 
are replaced by the A"*" projectors. The Lorentz invariant functions Fe{Q^) 



and Fm{Q^) denote the electric and magnetic formfactor, respectively, which 
depend on the photon momentum Q^. When the left hand side of this equation 
is calculated with the above loop diagrams in the Breit frame, i^e(Q^) and 
FmiQ"^) can be extracted by taking appropriate traces. 

In Fig. 2 the numerical results for the e.m. formfactors are shown. We ob- 
serve that an internal diquark structure (leading to a narrower vertex function 
in momentum space) seems to be necessary in our approach because otherwise 
the variation of the formfactors with around is only weak, the e.m. radii 
are then far too smalU. While the violation of charge conservation due to the 
impulse approximation can be compensated by choosing a suitable momentum 
distribution in the loop integrals, the absolute values of the magnetic moments 
are not large enough which means that 1^ diquarks are needed for a phe- 
nomenologically satisfying description. Since our model is covariant, timelike 
formfactors are also accessible, which can be seen in the figure: The formfactors 
are continous at — 0. 

5 Pion-nucleon formfactor 

When following Mandelstam's prescripton it is also possible to couple an ex- 
ternal pion current to the nucleon. Because of parity there is no coupling to 
the (scalar) diquark, so there is just one diagram to calculate (in analogy to 
equation (^TJ) but where the quark-photon vertex is replaced by the quark- 
pion vertex). The quark-pion vertex ( = pion Bethe-Salpeter amplitude) in 
the chiral limit (m^ — 0) is determined by the scalar quark selfenergy 

Tl{P^^Q,p')^ro^^r\ (24) 

tion of the pion-nucleon matrix element is given by 

J^{Q^) = {N{Pf,SfmN{P,,S,)) 

- u{Pf,Sf)[-i^g^NN{Q^)T'']u{P,,S,). (25) 

In analogy to the electromagnetic formfactors, g-nNN on the right hand side 
of this equation can be extracted when the left hand side is evaluated with 
the diagram where the pion couples to the quark and the diquark remains a 
spectator. Our numerical results for g-„NN are displayed in Fig. 3. Again a 
diquark structure is unavoidable to get the expected variation of the formfactor 
with the pion momentum. The absolute value of Qt^n n at the pion mass shell 
{Q^ = in the chiral limit) is in remarkable agreement with the experimental 
value of about 14. 
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Figure 3: Numerical Solutions for the Pion-Nucleon coupling g-nNNiQ^)- The thin curve 
denotes our results with the diquark formfactor, the thick curve our results without a form- 
factor. 

6 Conclusion 

In this talk we have presented a covariant bound state approach for the nucleon 
where quarks and diquarks interact through quark exchange. Using an effective 
parametrization of confinement we calculated within the needed approxima- 
tions not only the nucleon mass and wave function but also electromagnetic 
and strong formfactors. ft has been discussed that an internal diquark struc- 
ture, entering through a diquark formfactor in the BSE is preferable to obtain 
a nucleon wave function which has a realistic width in momentum space. Fur- 
thermore the impleni(;ntation of axialvector diquarks is necessary to improve 
the magnetic observables. 

This work is just the first step towards a full satisfying baryon model 
of the proposed diquark-quark bound state; approach. Nevertheless the results 
obtained so far are quite encouraging. To benefit from the advantages of our co- 
variant and confining approach we will also apply this model to electromagnetic 
and strong NN* transitions, virtual Compton scattering and A-production. 
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